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ABSTRACT 

We explore the dependence of weak lensing phenomena on the background cosmology. We first gener- 
alise the relation between P^(o;), the angular power spectrum of the distortion, and the power spectrum 
of density fluctuations to non-flat cosmologies. We then compute for various illustrative models. A 
useful cosmological discriminator is the growth of P^ with source redshift which is much stronger in low 
matter density models, and especially in A-dominated models. With even crude redshift information 
(say from broad band colours) it should be possible to constrain the cosmological world model. The am- 
plitude of P^ (lo) is also quite sensitive to the cosmology, but requires a reliable external normalisation 
for the mass fluctuations. If one normalises to galaxy clustering, with M/L fixed by small-scale galaxy 
dynamics, then low density models predict a much stronger distortion. If, however, one normalises to 
large-scale bulk-flows, the predicted distortion for sources at redshifts Z s ~ 1 — 3 is rather insensitive to 
the background cosmology. The signals predicted here can be detected at a very high level of significance 
with a photometric survey covering say 10 square degrees, but sparse sampling is needed to avoid large 
sampling variance and we discuss the factors influencing the design of an optimum survey. Turning to 
weak lensing by clusters we find that for high lens redshifts [Z\ ~ 1) the critical density is substantially 
reduced in A models, but that the ratio of the shear or convergence to the velocity dispersions or X-ray 
temperature of clusters is only very weakly dependent on the cosmology. 

Subject headings: cosmology: theory - gravitational lensing - large-scale structure 



1. INTRODUCTION 

Weak lensing is the distortion of the shapes and sizes 
(and hence fluxes) of distant galaxies from tidal deflection 
of light rays by structures along the line of sight. This pro- 
vides a powerful probe of mass fluctuations on a wide range 
of scales from galaxy haloes (Valdes et al, 1984; Brain- 
erd et al, 1996) through clusters (Tyson, Valdes & Wenk, 
1990; Bonnet et al, 1993; Bonnet et al, 1994; Dahle et 
al, 1994; Fahlman et al, 1994; Smail et al, 1994; Smail 
and Dickinson, 1995; Tyson and Fisher, 1995; Fort et al., 
1996 Squires et al., 1996a; Squires et al., 1996b; Dahle et 
al., 1996; Fort et al, 1996 ) to supercluster scales (Valdes 
et al, 1983; Mould et al, 1994, Villumsen, 1996a). As the 
distortion depends on the distance to the sources — being 
generally larger for more distant sources — weak lensing 
also provides a way to constrain the redshift distribution 
for very faint galaxies (Smail et al, 1995; Luppino and 
Kaiser, 1996; Fort et al, 1996). 

The first quantitative predictions for the distortion ef- 
fect were made by Gunn (1967) who was primarily inter- 
ested in what limits this placed on the classical cosmologi- 
cal tests. Dyer and Roeder (1974) estimated the distortion 
using "swiss cheese" models — with results quite similar 
to modern estimates (see below) — though they dismissed 
the effect as too small to have much observational signifi- 
cance. Webster (1985) computed the increase in ellipticity 
of distant objects in models with strong small scale mass 
inhomogeneity. This is related to the weak lensing effects 
discussed here, though in fact the broadening of the dis- 
tribution of cllipticities only appears at second and higher 
order in the gravitational potential $, and vanishes in the 



'weak lensing' limit considered here. In weak lensing, one 
considers the galaxy ellipticity to be a 2-component vector- 
like quantity e a = e{cos2<p, sin2<p}, where tp is the posi- 
tion angle. The expectation value of e a vanishes in the 
absence of lensing, and one searches for a coherent statisti- 
cal anisotropy of the e Q distribution caused by intervening 
matter. 

Quantitative predictions for weak lensing from modern 
ab initio models for large-scale structure were made by 
Blandford et al (1991) and Miralda-Escude (1991) who 
computed the rms distortion or 'polarisation' of distant 
galaxy shapes and found that it could be quite large — 
on the order of a few percent in rms shear — for the then 
popular theoretical models such as CDM, and that the 
rms shear increased as the 3/2 power of comoving source 
distance. These works also estimated the shear autocor- 
relation function. In Kaiser (1992, hereafter K92) this 
analysis was generalised to include a distribution of source 
distances, and extended in a number of other ways: It was 
shown how the angular power spectrum of the shear field 
was related to the power spectrum of the density fluctua- 
tions in three dimensions; that shear at observable levels 
on degree scales was an inevitable consequence of density 
fluctuations inferred from large-scale deviations from Hub- 
ble flow, or 'bulk-flows' (see Strauss and Willick, 1995 for 
a detailed review); that the (fourier transform of the) pro- 
jected surface density could be computed from the shear; 
and also that one could determine the M/L for foreground 
structure by cross correlating the foreground galaxy den- 
sity distribution with the shear. However, this analysis, 
like those of Blandford et al, and Miralda-Escude, was 
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restricted to the Einstein - de Sitter cosmological back- 
ground. More recently, Villumsen (1996c), and Bar-Kana 
(1996) have discussed the distortion in open cosmologies. 
Here we will explore in more detail how the cosmological 
background affects weak-lensing observables. In §|| we give 
a brief yet self-contained derivation of the distortion ten- 
sor — which describes the mapping between angles on the 
sky at the observer and position on some distant source 
plane — as a projection of the transverse components of 
the tidal field along the line of sight. In §|| we consider 
lensing by large-scale structure. We first derive the rela- 
tion between the angular power spectrum of the distortion 
and the power spectrum of the density fluctuations, and we 
then consider various illustrative models for P{k) and then 
discuss the feasibility of these observations and sampling 
strategy issues. In §^ we consider lensing by individual 
clusters. Our main goal is to elucidate the dependence of 
weak lensing phenomena on the background cosmologies; 
this will enable us to understand to what extent our con- 
clusions about the mass distribution and the redshifts of 
faint galaxies are cosmology dependent and, to turn the 
question around, to see to what extent weak lensing ob- 
servations, perhaps combined with other observations, can 
be used to constrain the cosmological world model. In a 
recent study, Bernardeau et al. (1996) have also consid- 
ered some aspects of weak lensing considered here, though 
their work emphasised more the possibility of measuring 
higher (than second) order moments. 



which, on scales much less than the curvature scale, 
and in cartesian coordinates, becomes ds 2 = a 2 (r] a f3 — 
2$8 a p)dr a drP where rj a p — diag{— 1, 1, 1, 1} which is the 
usual weak-field solution for a source Sp(r) where $ satis- 
fies Poisson's equation 



V 2 $ = 4irG5p 



(5) 



and where the Laplacian is taken with respect to proper 
distance. 

Photon trajectories in the spacetime are solutions of 
the geodesic equation 



2^a 



d 2 r 



dX 2 
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dr 1 " dr u 



9 ^~2 M )^XlX (6) 



To zeroth order in $, and for null geodesies, the time com- 
ponent of (|) is d 2 t/d\ 2 = ~H(dt/dX) 2 which one can 
solve to obtain the affine parameter dX — adt. Let us 
consider rays confined to a narrow cone around the polar 
axis: 9 <C 1, and let da 2 = dO 2 + d6 2 with 9 X — 9 cos <p, 
9 y = 9 sin ip. The zeroth order solution of the radial com- 
ponent of (||) is dz/drj = 1 where the conformal time is 
defined as usual by 77 = J dt/a. The angular components 
of (||), up to first order in $, d9ijdr\ are 



d 2 9 % 
dr/ 2 



= -2 



cosh z d9j 



9$ 



sinh z dr] sinh 2 z d9i 



(7) 



2. WEAK LENSING IN A FRW COSMOLOGY 

We take as the metric for the homogeneous and isotropic 
FRW background 

ds 2 = g a(i dr a dr S3 = -dt 2 + a 2 (t)(dz 2 + sinh 2 zda 2 ) (1) 

where da 2 = d9 2 + sin 2 9dip 2 . (Here z measures comoving 
separation; we will use uppercase Z to denote rcdshift). 
This is for the open case; for the closed world model we 
replace sinhz with sinz. The curvature radius a obeys 



H a Q 



a, 



ao 
a 



n A - 



± 1 



(2) 



or, in terms of the transverse comoving displacement of the 
ray from the polar axis measured in units of the curvature 
scale; X4 = ^sinhz, 



2d t $ 



(8) 



where di = d/dxi and dot denotes differentiation wrt con- 
formal lookback time z. The first term describes the ten- 
dency for neighbouring rays to diverge due to the hy- 
perbolic geometry, and becomes negligible in the limit 
Oo — * lj while the extra forcing term is, as usual, just 
twice the transverse gradient of the Newtonian potential. 
The general solution of m) is 



with positive/negative curvature term for open/closed 
models, and where H = (da/dt)/a is the expansion rate, 
a subscript '0' denotes the present value and where fi m , 
f^A are understood to be the present values of the density 
from matter and from the cosmological constant in units of 
the critical value p c = 3H 2 /(8irG). Evaluating (||) at the 
present we find that the curvature radius, which we shall 
use as the scale factor, is related to the current expansion 
rate by 

a o = 1^=^= (3) 

#ov 1 ~ ^0 

where fio = + ^A- We will focus on flat and open 
models (the former being considered as the limiting case 
of the latter), though it is straightforward to generalise the 
formulae below to the closed case. 

With the addition of small-scale matter inhomogeneity 
we can take the metric to be 



ds 2 



(1 + 2$)dt 2 + (1 - 2<&)a 2 {t){dz 2 + sinh 2 zda 2 ) (4) 



Ai sinh z + Bi cosh ; 



dz'<9 4 $(z')sinh(z-z') (9) 



which one can readily verify by direct differentiation. The 
constants of integration A and B are set by the boundary 
conditions. For a ray which reaches the observer (who we 
shall place at the origin of our coordinates) from direction 
9oi these conditions are Bi = 0, and Ai = 9 0i . 

If we consider a pair of neighbouring rays, and assume 
continuity of the potential <&, we obtain the geodesic de- 
viation equation 



Aa;, = Axi — 2Ax j d 1 d i ^ 



(10) 



This admits a solution as a perturbative expansion in $ 
where the n-th order term satisfies 



Ax, 



(«) 



Ax\ n) -2Axf- 1] d 3 d^ 



(11) 
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Thus, starting with the zeroth order solution Ax^ = 

A#jSinhz, one can obtain the solution for Ax^ in the 
form of (KJ), which can then be used to give the forcing 
term for the next approximation and so on. Here we will 
restrict attention to the first order solution: 



z 

Axi = A9i sinh z - 2A0j J dz' sinh(z') sinh(z - z')djdi$ 

o 

(12) 

This corresponds to evaluating the forcing term using the 
zeroth order separation; this being valid either in the 
'weak-lensing' approximation where the geodesic devia- 
tions are a small perturbation, or in the 'thin lens' ap- 
proximation, where the ray focusing may become large at 
great distances but where the change in the separation of 
the rays as they pass through the lens is small. 

Equation (Jl3) gives the mapping between angles at the 
observer and distance on some distant source plane at z s : 



Axi(z s ) = (Si m - ipi m )smhz s A9r, 
where we have defined the distortion tensor 



/ , sinh z smh(z s — z) „ „ 
sinh z s 



(13) 



(14) 



in agreement with Bar-Kana (1996). This is an observ- 
able quantity; the traceless parts of ip causing distortion 
of shapes of distant galaxies and the trace causing amplifi- 
cation and hence modulation of the counts of galaxies. In 
reality, we deal with the mean distortion tensor averaged 
over n(z), the distribution of distances to the galaxies 

iplm = J dz 8 n(z 8 )ipi m (z s ) = J dzg(z)did m $ (15) 

where 



g(z) = 2 sinh z I dz'niz') ^^- ,^ (16) 

sinh z 



is a bell-shaped function which peaks at roughly half of 
the background source distance, and where we have nor- 
malised n(z) so J dzn(z) = 1. 



3. LARGE-SCALE STRUCTURE 

We now consider lensing by large-scale structure. We 
first derive an expression for the angular power spectrum 
of some projected quantity (be it galaxy counts, image dis- 
tortion or whatever) and the corresponding spatial power 
spectrum. This is the fourier space analogue of Limber's 
equation (K92), here generalised to hyperbolic geometries. 
We then consider various illustrative models for P(k) of 
increasing degrees of realism and then discuss the fea- 
sibility of these observations, sampling strategy and the 
prospects for probing large-scale structure via the ampli- 
fication rather than shear. 



3.1. Limber's Equation in Fourier Space 

A common problem in astronomy is that one observes 
some quantity on the sky which is the projection of some 
three-dimensional random field or point process, and one 
would like to infer the statistical properties of the latter 
from the former. An example is galaxy clustering, where 
one would like to relate e.g. the angular correlation func- 
tion of the galaxy counts w g (9) to the spatial correlation 
function £ g (0). The solution to this was given by Limber 
(1954). Let the projected field be 

F(9) = J dzq(z)f(6 x sinh z,0 y sinh z,z) (17) 

where / is the spatial field written as a function of comov- 
ing coordinates (all in units of the curvature scale) and 
q(z) is some radial weighting function. The angular two 
point function of F is 

w F (A9) = (F(9)F(9 + A9)) = J dz J dz'q(z)q(z')(ff) 
~ J dzq 2 (z) J dz'£f(A9 x sinh z', A9 y sinh z' , z'\ z) 

(18) 

where £/(r; z) is the spatial two-point function of the field 
/ at lag r and conformal lookback time z, and we have as- 
sumed that q(z) is slowly varying compared to the scale of 
the density fluctuations of interest and also that these fluc- 
tuations occur on a scale much smaller than the curvature 
scale. This is Limber's (1954) equation, which expresses 
as an integral of the spatial two-point function, If we 
fourier transform ( p"8"| ) we obtain the angular power spec- 
trum Pp(ui). If we define the transforms 



F(uo) = f d 2 9F(9)e- luJ 
/(k)=/d 3 r/(r)e- lk -' 



(19) 



then under the assumption that the fields F, f are sta- 
tistically homogeneous (or more specifically that the two 
point function £f — (/(r)/(r')) depends only on separa- 
tion r' — r) we have 



<F(a,)F>')> 
</(k)/*(k')> 



(27T) 3 (5(k- 



- lo')P f {lu) 
k')P/(k) 



(20) 



where Pf{<->j) and P/(k) are the transforms of wf{9) and 
£/(r), so from (|lg| ) 

P P (w) - Jd 2 9w F (9)e— e = J ^P/tkJ/^W 

X J cp9e~ i ( UJ:c ~ k:c sinllz )^g _, '( IJ s _t » sinhz)0 B J ^1 g—ikzZ 1 

(21) 

The angular and z' integrals here are <5-functions 
which pick out the particular spatial frequency k = 
{u x / sinh z, ujy/ sinh z, 0} which contribute to the angular 
power at frequency to = {w x ,w 9 }, and then invoking the 
assumed statistical isotropy of P/(k) we have 



dz 



sinh 2 2 



P/(w/sinhz;z) (22) 



This is the fourier space version of Limber's equation, and 
is somewhat simpler than ( |l8| ) as it gives the angular power 
spectrum of F as a single integral of the spatial power 
spectrum of /, and provides the generalisation of (A9) of 
K92 to hyperbolic geometries. As in flat space, it can 
be thought of as a convolution in log-frequency space of 
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the three-dimensional power spectrum of /. Equation (22) 
can be used to relate the angular power spectrum of galaxy 
counts to the 3-dimensional spectrum of galaxy clustering, 
in which context F and / would be the density contrast 
of galaxies on the sky and in space, and q{z) would be the 
normalised distribution of galaxy distances n{z). Given 
a specific prediction for P(k), (g2|) enables one to predict 
P(w), or, given suffiently high signal to noise, one can 
deconvolve P(k) from P{lo) (see Baugh and Efstathiou 
(1994), who used this to extract the three dimensional 
power spectrum of galaxy clustering P g (k) from the angu- 
lar power spectrum P g (w) from the APM survey). 

What is the advantage of angular power spectrum analy- 
sis (PSA) over the angular auto-correlation function? One 
minor advantage, as we have seen, the former is somewhat 
easier to compute from P{k) 1 particularly when there is 
evolution of P(k). The real advantage of PSA for galaxy 
clustering, however, is that it is easy to compute the real 
uncertainty in the power estimates and that the error ma- 
trix for P(w) is nearly diagonal; aside from a readily calcu- 
lable short range correlation one scales 8u ~ 1/6, where 
8 is the dimension of the survey, estimates of P(w) at dif- 
ferent frequencies are statistically independent. Neither of 
these pleasant properties hold for correlation analysis. In 
weak lensing as we shall now see there is a further advan- 
tage in that the observable is a symmetric 2x2 tensor, and 
there are various correlation functions one can form, and 
making sense of the inter-relation between these is much 
simpler in terms of power spectra. 



3.2. Power Spectrum of the Distortion 

It is now very easy to compute the power spectrum of 
the distortion P^(uj) in terms of the 3-dimensional density 
field power spectrum P(k), or equivalently, in terms of 
P$(/c), the power spectrum for the potential fluctuations. 

The distortion tensor ( |l5|) can be written as the second 
derivative of a 'projected potential': 

i>i m {0) = didm%{0) (23) 

where di here and henceforth denotes d/d9i and where 



$ p (6>) = j dz 9 ^ $(6 x smhz,6 y smhz,z) 
J sinh z 



(24) 



which is in the form of dl7| ) with q = g/ sinh 2 z. In fourier 
space, differentiation is equivalent to multiplication by iuj: 
di^p — * iu>i$p so ipi m (uj) — ~u>iUJ m $ p (uj) and therefore 
the two point function for the distortion is 

(ipij(uj)ip* m (u!')) = (27r) 2 £(u/ - w)CjiU)jUJi(b m P^(Lo) (25) 

which nicely factorises into a pure angular term involving 
the unit wave vector u> — which, as we will discuss below, 
can be used as a test of the integrity of the data — and 
the 'distortion power spectrum': 



P^oj) = uj 4 



dz 



9 2 {z) 
sinh 6 z 



P$(u>/ sinhz; z) 



(26) 



which is a function only of |u>|. Equation ( p6[ ) expresses 
P^{u) as a convolution in log- frequency space of the power 
spectrum of potential fluctuations Pj>(fc), or equivalently 



to the power spectrum of density fluctuations Ps(k) 7 which 
is related to P$(fc), through Poisson's equation (see (35) 
below). Villumsen (1996c) has obtained an expression sim- 
ilar to (|2^), but finds a different angular dependence. 

3.3. Models for P*(w) 

We now compute P*(w) for various models for the 3- 
dimcnsional power spectrum P$(fc). We first consider the 
effect of power concentrated at a single spatial frequency. 
We next consider power law models and finally we con- 
sider empirical models which seem to fit most of the data 
on galaxy clustering as constructed by Peacock (1996). 

3.3.1. S-function P(k) 

To explore how the background cosmology affects the 
interpretation of P^(uj) let us first compute the angular 
power spectrum for a narrow band of power at some fre- 
quency with present physical wavenumber fc p h ys = k*, 
i.e. Pd>(k) — 27r 2 ($ 2 )(a fc») _2 (5(fc — a fc*). Inserting this 
is (|26|), and considering the effect on sources at a single 
rcdshift (n(z) — * S(z — z s ) so g(z) — ► 2sinhzsinh(z s — 
z)j sinh z s ) we find 



P^w) = 87r 2 ($ 2 )w 



/ 2 (z)sinh 2 (z s - z) 
sinh z, coshz 



(27) 



where z = smh — (u>/(aoA;*)) and where P^ vanishes for 
oj > aofc, sinhz s . Here / = (clq / a){5 / 5q) is the growth fac- 
tor for the potential, expressed as a function of lookback 
time. Equation (^) is calculated as follows: Equation (§) 
gives H/Ho as a function of l + Z = ao/a. Starting at some 
very high redshift we compute r\ — J dt/a = J daH/a to 
obtain the conformal lookback time z — r/o — r/ as a func- 
tion of Z. We also integrate the equation for the growth 
of density perturbations: 



dt 2 



2H- 



,d5 

~dt 



Gp m 5 = 



(28) 



to obtain S/Sq for the growing mode. The result is shown, 
for three representative cosmological models in figure ([!]). 
In all cases, the angular power spectrum is a bell-shaped 
curve, peaking at roughly half the maximum angular fre- 
quency w m ax = aofc*sinhz s . Also, the total power in- 
creases quite rapidly with increasing source redshift, this 
trend being strongest for the A-dominated models. 
In figure (0) we show the total power 



(29) 



/g^( W )=47r<* a >(&) (l-^o)- 3/2 



2 / \ sinh 2 z sinh 2 (z s —z) 
~~ sinh 2 z s 



x / dzf 2 {z) — r- 





and also the mean (power weighted) angular frequency: 



J dujiJ 1 P^,(uj 
J duiuiP^ui) 



(30) 



as a function of source redshift. At low source redshift 
the distortion power spectrum is essentially independent 
of the background cosmology, as one might expect, but at 
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high Z s the low matter density models predict higher dis- 
tortion. This is due in part to the fact that <J> decreases 
with time in these models, and in part to the greater path 
length back to a given source redshift. 
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Fig. 1. — Angular power spectra for a delta-function spatial power 
spectrum and for various source redshifts. We have used uiPj,(ui) /2ir 
as ordinate (with linear frequency as the abscissa) so that the area 
under the curves represents the total power. The solid line is the 
EdS model, the dashed line is an open Q m = 0.2 model and the 
dash-dot line is a A dominated model also with Q m = 0.2. The nu- 
merical values are for (4 >2 ) = 1, fc* = Hq, so for any other values the 
vertical and horizontal scales should be multiplied by (<3? 2 }(fc»/Ho) 2 
and k t /Ho respectively. 




Fig. 2. — The upper panel panel shows the total distortion power 
(for a 5-function potential power spectrum) as a function of the 
source redshift. As before, the numerical values assume ("l? 2 ) = 1, 
fc„ = Ho, and one should multiply the vertical scale by (<E> 2 )(fe*/_ffo) 3 
for other values. The lower panel shows the power weighted mean an- 
gular frequency for a fiducial physical wavenumber k„ = Ho, which 
shows that in low matter density models the power from a given 
physical scale appears at a somewhat (up to about a factor two) 
larger angular frequency. More striking, however, is the difference 
in the total power between the models; low matter density models, 
and especially A-models, predict much stronger distortion. This is 
assuming the same amplitude of 3-D potential fluctuations for the 
different cosmologies. The motivation for this and alternative nor- 
malisations are discussed in the text. 



3.3.2. Power-Law P{k) 

These results for a 8- function P$(fc) are most useful to 
show what spatial scales we are probing when we mea- 
sure the angular power at some frequency. For a realistic 
P<s>(k) we will see a blend of angular spectra of the form 
(pTj), As a next step towards realism we now model P$(fc) 
as a power law P$(fc) oc fc"~ 4 (following the usual conven- 
tion that the density power spectrum scales as k n ). To set 
the normalisation consistently (corresponding to a given 
metric fluctuation variance on a given physical scale), let 
us take 



k 3 Rj(k) 
2tt 2 



( fcphys 



= ($ 2 



(31) 



G 



so (3> 2 )„ is the contribution to the variance of the poten- 
tial per log interval of angular wavenumber at some fidu- 
cial (physical) wavenumber k*, and which gives, for the 
distortion power per log interval of wavenumber, 



w 2 P,/,H/(27r) = 4irC n (z s )($ 2 



,2+n 



(32) 



where the cosmology dependence is all hidden in the func- 
tion 



C n (2 s ) = (l-fi ) (n - 1)/2 dzf 2 (z) 



sinh 2 (z s — z) 
sinh™ z sinh 2 z s 



(33) 



Note that for n = — 2, which gives equal variance in ip per 
log interval of wave-number, (^) and ( |33"| ) are essentially 
identical to (p9|). The dependence of P^Joj) on source red- 
shift is shown in figure (0) for various values of the spectral 
index n. 
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Fig. 3. — Dependence of the power (per log interval of wavenum- 
ber) on source redshift for our three example cosmological models 
now assuming a power law spectrum of density fluctuations. The 
quantity plotted here is actually AnC n , and should be multiplied by 
{<S> 2 ),(k,/H ) 1 - n w 2+n to obtain the real power per log interval of 
angular wave-number. 



3.3.3. Normalisation 



It is readily apparent that the distortion power is a 
strongly increasing function of source redshift in all cos- 
mologies (though the trend is n-dependent). We also see 
that the predicted distortion for high-Z sources is sensi- 
tive to the cosmology, with low matter density models, and 
A-models in particular, giving a much stronger predicted 
signal. This might seem to be at odds with the conclusions 
of Villumsen (1996a) and Bernardeau et al. (1996), both 
of whom find the distortion to be a strongly increasing 
function of Q m (Bernardeau et al. (1996) find variance for 
the top-hat averaged shear approximately proportional to 
fi^ 5 for instance). However, the difference is simply one of 
normalisation; we have normalised to a given rms potential 
fluctuation whereas Villumsen and Bernardeau et al. have 
implicitly normalised to a given rms density contrast, and 
had we done this it would have reduced our predictions 
by a factor Ct^. The question of what is the appropriate 
normalisation is an interesting one, to which there is as yet 
no completely definitive answer. Were it the case that all 
one knew about the universe was the galaxy correlation 
function £ g , then one could reasonably make a case for 
normalising to a given density contrast; implicitly assum- 
ing, in the absence of any evidence to the contrary, that 
galaxies are unbiased tracers of the mass. However, there 
is actually a wealth of data from dynamical studies of var- 
ious kinds to suggest that the appropriate normalisation 
for the density contrast is quite strongly fi m dependent; 
in high Q models the galaxy distribution must be biased 
and it makes sense to fold this into ones predictions for 
the rms distortion. 

One line of evidence comes from small-scale pairwise 
velocities; in the regime where these motions are in equi- 
librium these measure the present potential fluctuations 
directly, and would therefore lead one to normalise as we 
have done. Thus, if one assumes a scale independent bias, 
with M/L tied to small scale 'cosmic virial theorem' mea- 
surements (Davis and Peebles, 1983) then, as shown in 
figure ||, low fl m models then predict much stronger dis- 
tortion at high Z s (this assumes that one is using 'real- 
space' estimates of £ g ; if one uses redshift survey based 
estimates then one should allow for the boosting of power 
due to streaming motions (Kaiser, 1987)). Another line of 
evidence comes from peculiar velocities on large scale; the 
so-called 'bulk flows' (e.g. Strauss and Willick, 1995 and 
references therein). The advantage of these observations 
is that they directly probe the mass, and give a normalisa- 
tion to the mass fluctuations on large-scales where, as we 
shall see, we expect weak lensing to be most powerful. The 
disadvantage is that they are very hard to measure reliably 
and estimates of the mass power spectrum derived there- 
from have very large 'sampling uncertainty' as the obser- 
vations typically probe only a few independent fluctuation 
volumes. Here the velocity on a given scale is, crudely 
speaking, given by the potential gradient times the age of 
the universe, and as the latter is greater in low Q m models 
we would infer a lower potential fluctuation amplitude for 
a given amplitude for the streaming motions. A detailed 
analysis shows that this reduces the predicted distortion 
by roughly a factor fJ^ 4 for low fl m models (in amplitude, 
that is, corresponding to a factor fi^ 8 in power). The 
predicted distortion power is shown in figure [| for this 
normalisation. Now we find that the greatest distinction 
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between models appears at very low source redshift, but 
for sources at Z s ~ 1 — 3, and realistic spectral indices, 
the predictions are only mildy dependent on cosmology 
and hence on O m . A third line of evidence comes from 
the abundance of clusters as a function of velocity disper- 
sion or X-ray temperature. With e.g. the Press-Schechter 
(1974) model for the mass function, these data can be used 
to normalise the amplitude of mass fluctuations (Cole and 
Kaiser, 1989). In high Q models, and with the best cur- 
rent data, this gives er 8 - 0.57f2~ a56 (White et al.\ 1993) 
which is very similar to the scaling from peculiar velocities; 
however, since low density models tend to produce power 
spectra with more large-scale power (as the horizon size 
at matter-radiation equality is increased), the predicted 
power at large scales in low-f2 m models would not be sup- 
pressed as much as with the bulk-flow normalisation and 
one would expect to find something intermediate between 
the predictions shown in figures |^, 
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Fig. 4. — Predicted distortion power versus source redshift as in fig- 
ure ^ but now normalised assuming a fixed amplitude of bulk-flows. 



There is still considerable uncertainty in the measure- 
ments mentioned above and in their interpretation, and so 
there is considerable slop in the predicted distortion. How- 
ever, we would stress that with any of the normalisation 
methods described above, and for sources at reasonably 
high redshift Z s ~ 1 — 3 say, and for plausible spectral 
indices n in the range —1 to —2 say, the A model pre- 
dictions are at least as high as in the Einstein - de Sitter 



model, and in open models the predicted distortion power 
is reduced by at most a factor two and certainly not by a 
factor 



3.3.4. Empirical Models for P{k) 

As a final example, we give the distortion power pre- 
dicted according to Peacock's (1996) model for the linear 
power spectrum, which, for high-f2 m incorporates a mild, 
though plausible level of bias, and gives a very good fit to 
most galaxy clustering data. Poisson's equation (H) gives 
P*(fc) = (47r) 2 a^- 4 P,(fc) = (9/4)fi*,*- 4 (l - Q )^Ps(k) 
where Ps(k) is the power spectrum for (mass) density con- 
trast, in terms of which (Eq) becomes 



2w 3 



/Cp hys P|' hys (fcphys)/27r 2 is the density contrast 



dz 



g 2 (z) sinh z/ 2 A 2 
(l-^o) 2 



(34) 



where A 2 = 

power per log interval of wavenumber and should be eval- 
uated at fcpnys — wPoVl — sinhz, and at conformal 
lookback time z. For sources at a single redshift this be- 
comes 



187r 2 fii 



dz 



sinh 3 z sinh 2 (z s — z)f 2 A 2 
(1 — fio) 2 sinh 2 z s 



(35) 



This, or more generally (|T^) and ([54j), is the most con- 
venient form if one wishes to predict the distortion from 
e.g. COBE normalised ab initio models. The results are 
shown in figure @ for Peacock's model for A 2 (fc), which are 
quite similar on form to MDM models. Again, we see that 
the difference between the different cosmological models, 
when normalised realistically, is quite mild. The quantity 
plotted here is A\, the contribution to the variance of the 
trace of ipi m per log interval, which is four times the vari- 
ance in the shear 7 or the convergence k, so these models 
are predicting rms shear, convergence ~ 1% for u> > 100. 
This is similar to the predictions of K92, though in fact 
the mass fluctuations assumed here are somewhat lower 
while the adopted redshift is higher to reflect the growing 
evidence (both from lensing and spectroscopy) for a signif- 
icant high-redshift component at the relevant magnitude 
limits. 

The formalism we have developed allows us to com- 
pute the corresponding quantities for any given distribu- 
tion function for the background galaxy distances, but, 
generally speaking, the results are very similar to the sin- 
gle source plane results for a single plane at the median 
distance (see e.g. K92 for examples). 
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Fig. 5. — Distortion power according to Peacock's model for P(k). 
As before, the solid, dash and dash-dot lines are for EdS, open and 
A models, but here the low density models have Q m = 0.3. These 
predictions are based on a linearised power spectrum. This should 
be valid at large-scales (the strongest distortion here derives from 
fluctuations with wavelengths ~ lOO/i Mpc and are quite accu- 
rately linear) but will tend to underestimate the distortion at small 
scales where non-linearity acts to boost the power considerably at 
late times. The straight solid line is the expected noise (1-sigma) 
due to measurement errors for a 3-degree square survey with realis- 
tic number density and intrinsic ellipticities as described in the text 
and for a resolution dlnuj = 0.25. The large error bars (arbitrarily 
attached to the open model prediction) illustrate the sampling noise 
for a 3-degree survey at this resolution. These can be reduced consid- 
erably by sparse sampling, as illustrated by the error bars attached 
to the A-model which are for a sparse sample of side 9 degrees. 



3.4. Feasibility and Strategy 

We have computed above the power spectrum for the 
distortion tensor tpi m . The quantities we actually measure 
are the convergence and shear 



where 



M llm = 



1 

-1 



M 2lm = 



1 

1 



(36) 



(37) 



The shear 7 a is measured from the shapes of galaxies, 
while the convergence k can be measured directly from 
the modulation of the counts. The latter tends to be rela- 
tively noisy (see Kaiser et at, 1994), so we will focus, for 
the moment, on the shear. 



3.4.1. Shear based P^(to) 

A procedure for estimating the power spectrum P^(uj) 
was outlined in K92. That analysis assumed a simple 
square survey geometry. Here we will generalise this to 
more complex survey shapes. Let us assume that we have 
observations of a set of N galaxies with positions {8 g } and 
that we have measured their shapes to obtain a set of prop- 
erly calibrated shear estimates, i.e. each galaxy provides 
an estimate ^ a = 7 a (6 g ) + 7^ nt where 7^ nt is a measure of 
the random intrinsic ellipticity of the galaxy plus measure- 
ment error. The first step is to take the fourier transform 
of this set of shear estimates: j a (u>) = la exp(— iuj ■ 6), 
which we can write as 



7a(«) = / d 2 9n{e) la (6)e 



-iuj-9 



(38) 



g ). Next we 



where we have introduced n(9) = Y 5{0 
convert this to an estimator of the dimcnsionless surface 
density R as follows: From (^6|) and ( p3|) it follows that the 
gradients of the surface density and shear are related by 
d[K = Maimdm'ja (Kaiser, 1995) and hence that V 2 k = 
Maimdidm^a, or, in fourier space k(lo) = M alm (biLb m ^ a (uj) 
which suggests the estimator 



k(uj) 



,7a(w) = C a (uj)%(uj) 



(39) 



where we have defined c a = {cos2ip, sin2t/?} and where tp 
is in turn defined by u>i — \ uj cos ip,uj sirup}. If one takes 
the inverse transform of ( |39| ) one obtains the surface den- 
sity estimator of Kaiser and Squires (1993). Here we will 
use (p9h t o obtain an estimate of the power spectrum P^. 
FmmlMj, (H), @ and © we find 



(|«H| 2 ) =iv( 7l 2 ) + CQ H C/3 H 

x / ^^c a (u/) C/3 (u/)|n(u; - u/)| 5 



(40) 



so (|k| 2 ) does indeed provide an estimate of the power con- 
volved with |n(w)| 2 plus a constant 'shot noise' term N^f) 
where (jf) is the mean square shear estimate (per com- 
ponent) due to random intrinsic shapes and measurement 
error. 

For a uniformly sampled square survey geome- 
try of side O the convolving kernel is |n(cj)| 2 = 
TV 2 sine 2 [ui x 6/2) sine 2 (u> y 6/2) which falls off rapidly for 
uj 3> 6 (and one could imagine introducing a weight func- 
tion tapering towards the edge of the box to apodize the 
kernel further). For realistic spectra, and for estimates of 
the power at uj ^> 1/0 the convolution integral in (40) 
is then dominated by frequencies very close to the target 
frequency (within Slu ~ 1/B), and one can then to a good 
approximation remove the relatively slowly varying factor 
PipCaCp from within the integral to obtain 



(|*M| 2 )=7V 



nP^ (lo) 



(41) 
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where we have invoked Parseval's theorem and where n = 
N/Q 2 is the mean number density of galaxies per stera- 
dian. To obtain our final estimate of the power we subtract 
the shot noise term to obtain P^(uj) = 4(|k| 2 — N{^ 2 ))/Nn 
and then average these estimates over a shell of frequencies 
in w-space with some width du> = Lodlnu. 

It is worth noting at this point that if we replace c a in 
(|39| ) by c' a = {sm2<^, — cos 2^}, or equivalently apply the 
90-degree rotation {71,72} — > {72, —71} to the shear esti- 
mates, then the estimated power should vanish (aside from 
statistical noise). This is a reflection of the fact that while 
a general shear field has 2 real degrees of freedom only 
one of them is excited by lensing. Using c a in (|3^) effec- 
tively projects out the active component, while c' a projects 
out the sterile component. This 'rotation test' provides a 
useful check on the integrity of the data, as most sources 
of spurious image polarisation would be expected to excite 
both components. For further discussion of this see Kaiser 
et ai, 1994 and Stebbins, 1996. 

To estimate the uncertainty in P^,{u>) we need to make 
some further assumptions about the density fluctuations. 
We shall assume that the real and imaginary parts of k(u>) 
approximate a pair of gaussian random fields. This is cer- 
tainly appropriate for density fluctuations arising from in- 
flation, but is also a very good approximation for many 
highly non-gaussian (in real space) processes. An exam- 
ple is a shot noise process where even for a quite mod- 
est number of shots (say 5 or so), the transform of the 
shots becomes quite accurately gaussian as a consequence 
of the central limit theorem (see Kaiser and Peacock, 1992, 
for illustrative examples). The gaussian (transform) ap- 
proximation is also valid for density fields which are non- 
gaussian due to either non-linear gravitational evolution 
or biasing. Under this assumption we can calculate the 
two point function for the 'raw' power (k 2 (ui)k 2 (to + Aw)) 
and hence obtain the uncertainty in the shell averaged 
power. For a filled square survey, the two point function of 
the power is, like |n(w)| 2 , a rather compact function with 
width 8u) ~ 1/6, so estimates of the power at wavenum- 
bers separated by more than 5u> are statistically indepen- 
dent, and computing the variance in the power averaged 
over some shell becomes essentially a counting exercise; 
one computes dN u , which is the effective number of inde- 
pendent modes, and then divides the mean power (signal 
+ shot noise) by y/dN^ (see Feldman et ai, 1994 for an ap- 
plication of this method to galaxy clustering in the QDOT 
redshift survey). 

For a simple filled square survey geometry the number 
of independent modes in each shell is just 



(7? 



2\l/2 



0.40 (as obtained from typical cluster lensing 



ttlo dlnu 
(6uj) 2 



(42) 



where Slo = 2tt/Q is the fundamental frequency. As- 
suming zero signal one obtains a statistical uncertainty in 
uj 2 P^(uj)/(2ir) due to random intrinsic ellipticities 



of 



(43) 



As noted in K92, for high-.Z sources the statistical uncer- 
tainty is very small compared to the expected signal. For 
example, with n~2x 10 5 galaxies per square degree, and 



studies for integrations of a few hours on a 4m class tele- 
scope) and for frequency resolution of say din a; = 1/4, or 
four bins of power per log interval of frequency, we obtain 
ct(A 2 ) ~ 1.15 x lO~ 6 (27nj/0) which is shown in figure || 
for a survey of side = 3 degrees. 



3.4.2. Sampling Strategy 

The measurement noise estimate ([l3]), is a tiny (~ 1%) 
fraction of the power for our fiducial 3-degree survey field 
(assuming Peacock's estimate of the power), so one would 
expect, in such a survey, to detect the power at something 
like the ~ 100-sigma level. This is very nice. However, it 
does not imply a similar precision in determining the true 
mean Universal power. For low spatial frequencies the 
uncertainty in the ensemble average power will be dom- 
inated by the fact that we only have a small number of 
independent modes. For a 3-degree survey, the fundamen- 
tal frequency is Slo ~ 120, and the fractional uncertainty 
in P^(ui) at say 25ui would be around 50%, rather than 1%. 
This 'sampling uncertainty' is shown as the error-bars in 
figure |5| and clearly dominates over the measurement noise. 

The situation here is very similar to that in galaxy clus- 
tering studies from redshift surveys where even though one 
might have a sample of many thousands of galaxies, the 
number of independent structures on the largest scales is 
quite small. In this situation one can measure the power 
in one's sample volume to extremely high precision, yet 
the value need not be representative of the ensemble aver- 
age power spectrum. Whether this 'sampling uncertainty' 
is relevant depends on what one wants to use the data 
for. One useful application of redshift surveys is to apply 
the cosmic virial theorem and obtain an estimate of the 
ft of matter clustered like galaxies. This essentially in- 
volves taking the ratio of the pairwise velocity dispersion 
to the galaxy clustering amplitude (Davis and Peebles, 
1983). Both of these statistics may have large sampling 
noise, but their ratio (under the assumption that there is 
a universal mass per galaxy) is not affected by this. Simi- 
larly, in weak lensing, one can perform a cross-correlation 
between the shear of the faintest galaxies and the surface 
number density of somewhat brighter galaxies (chosen so 
that their n(z) peaks roughly half way to the background 
galaxies), to obtain an estimate of the mass to light ratio 
(K92). As with virial analysis, this does not involve the 
sample variance, and M/L can be measured to extremely 
high precision using a filled survey. If, however, one's goal 
is to use P^(lo) to distinguish between e.g. the three mod- 
els shown in figure then the sampling noise is clearly a 
serious handicap. 

As with galaxy clustering, sparse sampling (Kaiser, 
1986) could be quite helpful here. Imagine one were to 
observe a similar number of galaxies, but with Nf sparsely 
spaced fields of size / scattered over a much larger square 
(of side denoted by as before). The function |n(w)| 2 , 
which should be thought of as a 'instrumental point spread 
function' through which we measure the power will now be 
much narrower (by a factor ~ y/J, where / is the areal fill- 
ing factor of ones survey). This has two benefits; first, the 
fundamental frequency would decrease, so one would be 
able to probe beyond the peak in the spectrum (clearly an 
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interesting and cosmology dependent attribute). Second, 
in each region of frequency space one will obtain 1/ / times 
as many independent estimates of the power, so the sam- 
pling uncertainty decreases by a factor l/yj, which can 
be a considerable gain. 

There is a price, however, for this increased resolution 
and precision, which is 'aliasing'. For a sparse survey, the 
convolving kernel |n(w)| 2 will now have side-lobes in ad- 
dition to the central peak which extend to relatively high 
frequencies uit ~ 2tt/0/, and our power estimator |k(w)| 2 
will in general contain some contribution from lo' =/= lu. 
The structure of the side-lobes depends on how one lays 
out ones fields. If this is done in a random or semi-random 
manner (perhaps by choosing a set of fields which avoid 
few bright foreground objects) then the side lobes look 
like (the square of) a random gaussian field smoothed on 
a scale 8ui and the mean strength of the side-lobes is sup- 
pressed by a factor ~ 1/Nf as compared to the central 
peak. If the fields are laid out on a grid then the side- 
lobes are as strong as the central lobe, but are spaced on a 
grid of spacing y/NfSui and cover only a fraction / of the 
frequency plane. Some examples are shown in figure (||). 

There are two aspects to the 'aliasing problem'; the first 
is mixing of power from frequencies similar to the target 
frequency: u' ~ u. For a random or grid-like survey this is 
small provided Nf 3> oj 2 Q 2 , and this condition says that 
one must sample several fields per wavelength of interest. 
If the fields are laid out in a line, a substantial mixing 
of power is unavoidable (as is the case for pencil beam 
or 2- dimensional redshift surveys). This is not absolutely 
disastrous, as one can always convolve one's theoretical 
predictions, but it seems an unwanted and unnecessary 
complication, and we would discourage this. Also, in order 
to apply the 'rotation test' described above it is necessary 
one has proper 2-dimensional sampling and that the above 
condition be satisfied. 



o 




Fig. 6. — The panels on the left illustrate three possible ways one 
might lay down fields for a weak lensing survey. In each case there 
are Nf = 36 fields, and, assuming a field size of 0.5 degrees the box 
side would be about 10 degrees. The top panel shows a filled sur- 
vey (3 degrees on a side), the middle panel shows a sparse survey 
with randomly placed fields while the bottom panel shows a 'VLA' 
style layout. The panels on the right show the corresponding 'point 
spread functions' for power spectrum estimation. The contours are 
logarithmic at 2 _n times the peak value. The filled survey has a 
wide central lobe, and hence very poor precision as the sampling of 
the power spectrum is very coarse (as discussed in the text, the vari- 
ance in the power is inversely proportional to the area of the central 
lobe), but the aliasing of high frequencies is very small. The random 
sparse survey has a very tight central lobe, and would yield a factor 
~ 3 improvement in precision over the filled survey. There is also 
rather good rejection of neighbouring frequencies close to the target 
frequency, but one can see the sidelobes (the typical height of which 
is suppressed by a factor 1/Nf relative to the central lobe) which 
extend to frequencies of order the inverse field size. For the models 
discussed here the power aliased through these sidelobes would be 
very small. The 'VLA' style survey has very poor behaviour and 
combines the worst features of the two strategies above; the cen- 
tral lobe is not much narrower than the filled case, which yields 
large sampling variance, while the side lobes are at least as strong 
as for the random sparse survey. Moreover, there is poor rejection 
of frequencies close to the target frequency, so the measured power 
spectrum at low frequencies will be quite distorted and one would 
need to apply deconvolution. 



The second aspect is aliasing of power from higher fre- 
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quencies ui' 3> u>. This is somewhat more model depen- 
dent as it depends on how the power varies with frequency. 
One finds from (40) that the aliased power will be small 
compared to the power one is measuring from the tar- 
get frequency, provided N f > w 2 9 2 A^(u;')/A^(a;). This 
again is physically reasonable; Nj /(cj 2 9 2 ) is the number of 
fields per square target wavelength, so when this condition 
is only marginally satisfied the small scale power induces 
a 'root-iV' contribution in the low-u k estimate which is 
equal to the intrinsic low frequency n fluctuations. Thus, 
if A 2 (ui) increases strongly towards small scales then one 
would want to increase the sampling density. However, if 
the models shown in figure || are a realistic guide A 2 ^ is 
actually falling with increasing frequency, so aliasing from 
small scales would not be a serious problem, and a rather 
sparse sampling would pay great dividends with little or 
no cost. It may be that this is overly optimistic, since, 
as mentioned above, non-linearity may boost the power at 
small scales at late times and cause the predictions shown 
in figure ^| to underestimate the lensing power. Nonethe- 
less, it should still possible to choose Nf so as to make the 
aliased power small. Assuming these conditions have been 
met, the power estimate is 



49 2 

P^) = _(|k| 2 -7V( 7i 2 )) 



(44) 



and the uncertainty in the distortion power per log interval 
of wave number is given by 



(if there is insufficient rejection of high frequency fluctu- 
ations one would need to estimate and then subtract an 
extra constant term from (fl4|), and there would also be 
an additional contribution to ( [45| ) from the fluctuations 
in the aliased pwer.) The two terms in ( (45| ) have op- 
posite dependence on the overall survey size 9, and the 
optimum choice (assuming a fixed instrument and inte- 
gration time) is when these two contributions are equal: 
9 op t = y/NP^p /4(7 2 ). For high redshift sources, and for 
large N, this would dictate a very sparse sampling indeed: 
/ = 4(7 2 ) / (nPip) ~ 10~ 2 for n ~ 2 x 10 5 per square degree 
and Z s ~ 3, but if one also wants to use the same survey 
to measure the shear for brighter and nearer objects one 
should adopt a correspondingly less sparse strategy. As a 
specific example, if one were to spread ones fields over a 
square say 3 times wider than the filled survey then the 
sampling uncertainty would decrease by a factor 3 while 
the 'shot-noise' term would increase by a factor 3. For 
Z s ~ 3 the former would still dominate, at least on large 
angular scales, and the resulting uncertainty is shown in 
figure U as the error-bars attached to the EdS spectrum. 
As one can see, the sparse sample has considerable im- 
proved precision and one also has better coverage of the 
behaviour of P^(cj) around the peak as the fundamental 
frequency is now reduced by a factor 3 also. To design 
a truly optimal survey for large-scale structure one really 
needs to know the true level of small scale power, so it 
would be prudent to first perform a smaller scale filled 
survey to establish this empirically. 



3.4.3. LSS and Amplification 



To close this discussion of large-scale structure, we con- 
sider what extra information can be gleaned from esti- 
mates of the surface density derived from the amplifica- 
tion. In general, we expect the rms convergence to be 
equal to the rms shear, and the convergence can be mea- 
sured directly as it causes a modulation of the counts of 
faint galaxies An/n = —2(1 — o)k, where a = dlnN(> 
l)/d\nl — 2.5d\ogN/dm&g is the logarithmic slope of the 
counts and the factor —2(1 — a) is known as the 'ampli- 
fication bias factor'. At the magnitudes relevant for weak 
lensing observations d log N/dm&g ~ 0.3 (we will discuss 
the dependence of a on waveband and colour presently) so 
large-scale structure will give rise to apparent clustering of 
galaxies on the sky with amplitude An/n ~ —0.5k. 

Now for a power-law spectrum of mass fluctuations in 
an Einstein - de Sitter cosmology the variance in k grows 
as z 1_n , where z is the mean comoving distance to the 
background galaxies (here z = 1 — (1 + Z) -1 / 2 ). On the 
other hand, if galaxies are unbiased tracers of the mass, 
then w(9) from the intrinsic spatial clustering decreases as 
(1 - z) 4 z~( 3+ ") (K92), so, in this model, the relative im- 
portance of the lensing induced clustering grows rapidly 
with increasing source redshift, and there is a crossover 
at Z s ~ 2 beyond which the induced clustering comes to 
dominate (K92, Villumsen, 1996a). As discussed above 
though, a mass-traces- light normalisation seems untenable 
for a high Q m universe and in cither a low Q m universe, or 
a biased f2 m = 1 universe, the expected induced clustering 
is sub-dominant. In any case, it is not entirely clear how 
one can separate the two effects observationally (without 
accurate redshift information). 

Another possibility is to cross correlate the surface den- 
sity k derived from lensing of faint galaxies with the surface 
brightness for a sample of brighter objects (chosen so that 
their n(z) peaks where the g(z) for the fainter background 
galaxies peaks). The procedure for doing this with k de- 
rived from the shear was described in K92, but this can 
equally well be done with n derived from the amplification. 
There is a slight complication here in that the predicted 
effect is an anti-correlation of background and foreground 
galaxy counts (as the bias factor is negative) , and this will 
be diluted (and perhaps overwhelmed) by contamination 
of the 'background' sample by faint foreground galaxies, so 
an accurate knowledge of the luminosity function is needed 
to properly account for this. 

For realistically normalised models we find rms k ~ 
1% on degree scales (for Z s ~ 3). For a shear based 
surface density estimate k 7 the uncertainty is <t(k 7 ) ~ 

V (7i) /N — 0.4/ViV and is much less than the expected 
signal for a survey like that discussed above. The un- 
certainty in the amplification derived surface density is 
a(K A ) = 0.5(1 - a)- VU + nP g )/N ~ 2y/{l + nP g )/N 
where P g is the power spectrum of galaxy clustering, which 
can be the angular power spectrum, in which case n should 
be the density of sources on the sky, or P g can be the spa- 
tial power spectrum, in which case n should be the three- 
dimensional number density, the result being the same. 
The ratio of the statistical uncertainty for the two tech- 
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niques is 



ct(ka) 
ct(k 7 ) 




(46) 



The factor 1 + nP g here (sometimes referred to as ~ 
(1 + AimJz)) is an effective clustering multiplicity. For 
a distribution of randomly placed clusters, it is just the 
number of galaxies per cluster, and quite generally, the 
factor nP g tells us by how much the variance of counts of 
galaxies exceeds that for a poisson distribution. 

For a power law angular correlation function w(9) = 
wo(9 1 '9o)~°' s , as seems to be a reasonable fit to the data, 
nP ~ 0.8(fc6» )~ 1 - 2 27rn'u;o6'o- There are now a number of 
empirical estimates of w(0) for faint galaxies which we can 
use to estimate 2irnwoOQ and some recent estimates (with 
#o = 1') are: (1.52, Couch, Jurcevic and Boyle, 1993; 1.1 
(B), 1.6 (R), Roche et ai, 1993; 1.45, Pritchet and In- 
fante, 1992; 1.74, Efstathiou et al, 1991; 3.8, Villumsen, 
Freudling and da Costa, 1996; 1.59, Brainerd, Smail and 
Mould, 1995). While these estimates span quite a range of 
magnitude limits (the number density varying from ~ 10 
per square arcmin for the brighter surveys to ~ 400 for 
the Hubble Deep Field) the estimates of 27rnwo^o are very 
stable and indicate nP ~ 2(wl') -1 - 2 . These results are 
discouraging in the extreme. Even on sub-arcmin scales 
where the clustering becomes negligible, the noise in the 
surface density inferred from amplification is already about 
5 times that for the shear based k estimates (so the ex- 
tra information contained in the amplification is meagre), 
and on larger scales the situation rapidly deteriorates. For 
cluster lensing (where one is probing structure on a scale 
~ 5' — 10') and where it might have been hoped that am- 
plification might resolve the 'mass-sheet degeneracy' prob- 
lem, clustering has already inflated the uncertainty by a 
factor 2 or so, so a cluster like A1689, which can be de- 
tected at the 10-sigma level in the shear, would only be de- 
tectable at the ~ la level (i.e. the diminution of the counts 
from lensing is about equal to the rms fluctuation expected 
from galaxy clustering alone). Worse still, for the degree 
scales of interest here u> ~ lOOrad" 1 ~ 0.03arcmin _1 and 
we expect nP to have grown to around 120, so the cluster- 
ing fluctuations exceed the poisson fluctuations by about 
an order of magnitude. We reach a very similar conclu- 
sion if we use Peacock's linearised estimate of P; which is 
not surprising as this power spectrum model is designed to 
fit the empirical data. Thus, on degree scales, one would 
expect ct(ka) ~ 50ct(k 7 ); so the extra information in the 
amplification based k estimate is negligible and, for sur- 
veys of the scale envisaged here, one would expect at best 
a marginal detection of the effect. The relatively high pre- 
cision allowed by the shear based surface density estimate 
relies on the assumption that the intrinsic shapes of galax- 
ies are uncorrelated, and could be compromised if in fact 
there are strong intrinsic alignments of galaxies on super- 
cluster scales. There have been a number of attempts to 
detect correlated orientations of galaxies in superclusters 
— with the hope of distinguishing between 'top-down' and 
'bottom-up' structure formation scenarios — but no con- 
vincing positive detections have been obtained (see Djor- 
govski, 1986 for a review), and there are no indications 
from weak lensing observations (from e.g. the rotation test) 
for any intrinsic alignments at problematic levels. 



The counts slope is somewhat waveband dependent (the 
counts in / are slightly flatter, and those in V somewhat 
steeper than we have assumed). This might suggest that 
one should use the redder passband, and it has been sug- 
gested (Broadhurst, 1996) that one would do better by 
selecting only red galaxies as they have an even shallower 
counts slope and hence a larger (negative) amplification 
bias. However, this does not help here as it now appears 
that faint red galaxies lie at relatively low redshift com- 
pared to their bluer cousins (Luppino and Kaiser, 1996) 
and this outweighs the gain from the flatter slope. It is 
somewhat unfortunate, though not entirely coincidental, 
that the most distant objects (the faint and blue galaxies) 
which suffer the greatest amplification have a very small 
amplification bias. 



4. CLUSTER LENSING 

We now consider lensing by an individual cluster which 
is assumed to dominate over the effect of foreground and 
background clutter. From (JlJ), and specialising to deflec- 
tions occurring in a single plane, we can write 



(47) 



where E is the surface mass density and E cr ;t is the critical 
surface density: 



iWl - ^o(l + Zi) sinhz, 
47rGsinhzi sinh(z s — z{) 



(48) 



which we have plotted against source redshift for a variety 
of lens redshifts in figure |7|. 

For low Zi the impact of cosmology is very weak, re- 
gardless of source redshift, but for Z s ~ 1 or higher we 
find a significant reduction in E cr ;t in the A dominated 
model (due to the increased path length) , and an observer 
living in such a universe but using the EdS E crit would 
overestimate the mass of a cluster at Z ~ 1 by almost a 
factor 2. 
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Fig. 7. — Critical surface density versus source redshift for various 
lens redshifts and for our three illustrative cosmological models. 



Fig. 8. — Relative distortion strength parameter /3 versus source 
redshift for various lens redshifts and for various values of f!o- 



It is also interesting to compare the mass inferred from 
lensing with that obtained from virial analysis and/or X- 
ray temperature information. Virial analysis gives 



£ = 



GOa(zi) sinh zi 



aH <7 2 VT~^(l + Zi) 
G sinh zi 



(49) 



where a is some number of order unity which accounts 
for the radial profile of the cluster; velocity dispersion 
anisotropy; departures from sphericity; departures from 
equilibrium; substructure; mass/light segregation etc. etc., 
so assuming that this can be done to sufficient accuracy, 
we should find 



nO sinh(z s - z t ) 

0{*l>*s)= „ ; „u_ ( 50 ) 



Airaa 2 



sinh z s 



This dimcnsionless ratio is an observable, and is depen- 
dent on the cosmology. To see how useful this is we have 
plotted (3 as a function of Zj, Z s for the three fiducial cos- 
mological models. It is clear that the effect of varying the 
cosmological parameters on this quantity is very weak. 
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5. SUMMARY 

The main new analytic result of this paper is the expres- 
sion (^6|) which gives the angular power spectrum of the 
distortion in terms of the 3-dimensional power spectrum 
of potential fluctuations (a somewhat more useful form 
for predicting P^(w) from models for the density contrast 
power spectrum P$ is given in (j35|)). We have used this to 
compute the distortion power for a number of illustrative 
models. 

We have shown that for a sharply peaked 3-dimensional 
spectrum, the scale on which the distortion appears is cos- 
mology dependent, being up to a factor ~ 2 smaller in 
low density models. This is reflected in the model predic- 
tions shown in figure |5| based on Peacock's models for P(k) 
which has a knee-like feature at A ~ 100/i _1 Mpc, and con- 
sequently the location of the corresponding peak in P^(lo) 
is comsology dependent. A more powerful cosmological 
discriminator is the growth of the distortion with redshift, 
this being stronger in low matter density models in general 
and in A-dominated models in particular. This requires 
that one have at least approximate estimates of the red- 
shift of the faint galaxies as a function of flux, but this 
should be feasible using approximate redshift estimates by 
fitting multicolour photometry to template spectra; Loh 
and Spillar, 1986; Conolly, et at, 1995; Sawicki et al, 
1996) This cosmological test does not require any external 
normalisation of the power spectrum, but does require that 
we should be able to measure the distortion with sufficient 
precision at both high and low redshift. 

One can also ask whether one can hope to pin down the 
cosmology by making use of external normalisation. This 
is certainly possible in principle, but it is not yet entirely 
clear what is the appropriate normalisation. Bernardeau 
et al. normalised to a fixed amplitude for the density con- 
trast and consequently found a very strong dependence of 
the predicted shear on the matter density: P^ cx ill n 5 or 
thereabouts. We have argued that this normalisation is 
unrealistic: if instead we normalise to galaxy clustering 
with a scale invariant bias and mass to light ratio fixed by 
small scale cosmic virial theorem analysis then we reach 
the opposite conclusion: low matter density models in gen- 
eral and in A-dominated models in particular then predict 
much stronger distortion at high redshift (the distortion 
being cosmology independent for low redshift). If on the 
other hand, one were to normalise to some value for the 
amplitude of large-scale bulk flows (still unfortunately a 
rather uncertain quantity) then oc £l"j 8 for very low 



source redshift, but for Z s ~ 1 — 3 and for realistic spec- 
tral indices n around —1 to —2 the predicted distortion 
is only very weakly cosmology dependent. A very similar 
result is obtained if one normalises to cluster abundances. 
This weak dependence on cosmology was also apparent 
when we computed the distortion power for Peacock's fit 
to galaxy clustering data, where all three illustrative mod- 
els agree in distortion power to within a factor 2. The high 
density model assumed a rather mild bias b — 1.6, and for 
stronger bias the difference between the models would be 
even less. 

With a realistic normalisation we predict rms shear at 
the ~ 1% level at degree scales for sources at Z s ~ 3, 
which should be detectable at the ~ 100-sigma level with 
a survey covering ~ 10 square degrees (which would con- 
tain ~ 2 x 10 6 galaxies). We found, however, that for a 
filled survey of this size the sampling uncertainty would 
much larger than the measurement noise, particularly at 
the largest scales which in some ways are the most inter- 
esting. For some applications the sampling uncertainty is 
irrelevant, but for the tests described above it is a seri- 
ous handicap. The sampling noise can be reduced con- 
siderably by adopting a sparse sampling strategy (at some 
small cost in increased measurement noise). An important 
constraint on the design of such sparse surveys is aliasing 
of power from small scales. While in principle this can be 
measured and subtracted to obtain a fair estimate of the 
true large-scale power, it is still an unwanted complication 
and, if the aliased power is dominant, the precision will be 
compromised. Peacock's empirically based models for the 
linear power spectrum predict very low power at high fre- 
quencies and would favour very sparse sampling (for deep 
surveys at least), but safest approach is to measure the 
high frequency power directly with a filled survey, and use 
this to determine the optimal sampling rate. 

Finally, we considered the impact of cosmology on mass 
estimates for individual clusters. We found that the crit- 
ical surface density was very similar in the matter domi- 
nated models, but is considerably lower for high redshift 
lenses in a A-dominated model. This is relevant to the 
result of Luppino and Kaiser (1996) who found a strong 
shear signal for the cluster msl054 at Zi = 0.83. In a 
A dominated model, the mass for this cluster would be 
reduced by about 40%. We also explored how the ratio 
of the lensing virial mass (or mass inferred from X-rays) 
depends on cosmology, but found this to be a very weak 
effect. 
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